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Abstract 

We introduce a new multipole formalism for polarized radiative transfer in 
general spacetime geometries. The polarization tensor is expanded in terms 
of coordinate- independent, projected symmetric trace-free (PSTF) tensor- 
valued multipoles. The PSTF representation allows us to discuss easily the 
observer dependence of the multipoles of the polarization, and to formulate 
the exact dynamics of the radiation in convenient 1+3 covariant form. For the 
case of an almost-Friedmann-Robertson- Walker (FRW) cosmological model 
we recast the Boltzmann equation for the polarization in to a hierarchy of 
multipole equations. This allows us to give a rigorous treatment of the gener- 
ation and propagation of the polarization of the cosmic microwave background 
in almost-FRW models (with open, closed or flat geometries) without recourse 
to any harmonic decomposition of the perturbations. We also show how ex- 
panding the intensity and polarization multipoles in derivatives of harmonic 
functions gives a streamlined derivation of the mode-expanded multipole hier- 
archies. Integral solutions to these hierarchies are provided, and the relation 
of our formalism to others in the literature is discussed. 

I. INTRODUCTION 

With the growing body of data relating to the anisotropies of the cosmic microwave 
background radiation (CMB), and the potential impact on cosmology of conclusions drawn 
from the analysis of CMB data, there is a strong case for developing a flexible, but physically 
transparent, formalism for describing the propagation of the CMB in general cosmological 
models. 

In a series of papers p|-[10[l, a 1+3 covariant and gauge-invariant formalism has been 



developed with a view to giving a model-independent framework in which to study the 
physics of the CMB. The approach is based on the projected symmetric trace-free (PSTF) 
representation of relativistic kinetic theory due to Ellis, Treciokas, and Matravers fill, and 
Thorne [0, and builds on the covariant and gauge- invariant approach to perturbations in 
cosmology (e.g. Ref. [0). Some of the benefits of the covariant approach, which were 
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emphasised in Refs. include: (i) clarity in the definition of the variables employed; (ii) 
covariant and gauge-invariant perturbation theory around a variety of background models; 
(iii) provision of a sound basis for studying non-linear effects; (iv) freedom to employ any 
coordinate system or tetrad. In its current form, the covariant formalism does not allow for 
the inclusion of polarization. It is this omission that we address here, by introducing a new 
multipole formalism which allows one to analyse CMB polarization in arbitrary cosmological 
models in the 1+3 covariant approach. 

We set up the multipole decomposition of the polarization and the exact equations of 
radiative transfer (the Boltzmann equation) in convenient 1-1-3 covariant form. The polar- 
ization multipoles are covariantly-defined PSTF tensors which leads to significant simplifi- 
cations in the algebraic structure of the scattering terms in the Boltzmann equation, and 
allows us to discuss easily how the multipoles transform under changes of reference frame. 
With the Boltzmann equation in 1+3 covariant form, the physical processes responsible for 
the production and evolution of anisotropics and polarization are particularly transparent. 
The PSTF multipole decomposition of the polarization suggests recasting the Boltzmann 
equation as a multipole hierarchy, in a similar manner to the treatment of the intensity in 
Refs. [0,0. This turns out to be quite involved, and we shall only present the results for 
the special case of an almost-FRW model here. (The exact, non-linear multipole equations 
will be given in a subsequent paper.) These multipole equations allow us to give a rigorous 
treatment of the generation and propagation of polarization in the CMB in almost-FRW 
models for all spatial geometries. 

Although the calculation of CMB anisotropics in an almost-FRW universe involves only 
linear perturbation theory, the complexity of the subject is increased when one allows for non- 



fiat spatial geometries, vector and tensor perturbations, and the inclusion of polarization [14 



To our knowledge, the only published formalism for handling the most general almost- 



FRW models is the total angular momentum method of Hu et al. [|I4[. In their approach, 
considerable simplifications result from the introduction of a normal mode expansion for the 
radiation where the local angular and spatial distribution is explicit. In the 1+3 covariant 
approach the local angular distribution of the radiation is analysed in terms of the PSTF 
tensor-valued multipoles, which can be done without recourse to a decomposition of the 
spatial distribution in harmonic functions. The multipole equations can then be used to 
analyse the evolution of anisotropics and polarization for a quite general perturbation. If 
required, the spatial dependencies can be handled by expanding the radiation multipoles 
in derivatives of harmonic functions, which gives a streamlined derivation of the mode- 
expanded multipole hierarchies. Combining the angular and spatial expansions we find a 
normal mode representation of the radiation which is equivalent to that of Wilson []T3[ for 
scalar temperature perturbations, while for vector and tensor modes we obtain the obvious 
generalisation of Wilson's method. 

The paper is arranged as follows. In Sec. we introduce the PSTF multipole decompo- 
sition of a polarized radiation field, relative to some timelike velocity field m", obtaining the 
1+3 covariant forms of the electric and magnetic parts of the linear polarization tensor. We 
also give an exact discussion of the transformation properties of the polarization tensor un- 
der changes of the velocity field m", and give non-linear expressions to first-order in relative 
velocities for the transformations of the polarization multipoles. In Sec. |II| we discuss the 
Boltzmann equation for polarized radiative transport. We give new, exact expressions for 
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the scattering term in the Thomson hmit and the integral solution for the linear polariza- 
tion, valid in general cosmological models. Specialising to almost-FRW models, we present 
the multipole form of the Boltzmann equation which shows explicitly that the evolution of 
the electric and magnetic parts of the linear polarization tensor are coupled through curl 
terms, much like the coupling of electric and magnetic fields in Maxwell's equations. In 
Sec. 1^ we decompose the multipole equations into harmonic modes for scalar perturbations 
and provide the integral solution which is central to the line of sight algorithm for efficient 
numerical evaluation of the CMB power spectrum [|1^. The analysis is repeated for tensor 
perturbations in Sec. [V], where a very direct derivation of the integral solutions is also given. 
In Sec. |V| we discuss the relation of our approach to others in the literature. Finally we 



close with our conclusions in Sec. [VI I| . An appendix summarises the PSTF representation 
of the scalar and tensor harmonics. 

We employ a (H ) signature for the spacetime metric Qab- Early lower case Roman 

indices a, b etc. refer to a general basis, with i and j referring to the 1, 2 components 
in an orthonormal tetrad, (cj)'*, i = 1,2. Round brackets denote symmetrisation on the 
enclosed indices, square brackets antisymmetrisation, and angle brackets the PSTF part: 
S(^ab) = {h^{J^h) ~ \hahh^'^)Scd, where the projection operator /i"* = g"-^ — u"-u^ with the 
fundamental velocity of the 1+3 covariant approach. The index notation Ai denotes the 
index string ai . . . a/, and the notation e^j denotes the tensor product e^^ . . . Cqj. We use 
units with c = G = 1 throughout. 



II. MULTIPOLE DECOMPOSITION OF THE RADIATION FIELD 

We describe observations from the viewpoint of an observer comoving with the funda- 
mental velocity field u°'. In a general cosmological model there is some freedom in the choice 
of although we must ensure that u"" is defined in a physical manner (such as by the time- 
like eigenvector of the matter stress-energy tensor, or the 4- velocity of some particle species) 
so that in the FRW limit, u'^ correctly reduces to the fundamental velocity of the FRW 
model. This restriction on u"- is necessary to ensure gauge- invariance of the 1+3 covariant 
perturbation theory. A photon with 4-momentum has an energy E and propagation 
direction relative to m'*, where 

= Eiu'' + e''). (2.1) 

For a given propagation direction e", the observer can introduce a pair of orthogonal polar- 
ization vectors (ei)" and {€2)°" which are perpendicular to and e": 

ntie^f = (el)^ (2.2) 

with a similar result for (62)". Here Tiab is the screen projection tensor that projects per- 
pendicular to both and e": 

'Hab = hab + CaCb. (2.3) 

We shall refer to tensors like (ei)", which are perpendicular to e"" and m", as being trans- 
verse. The set of vectors (ei)°, (62)", e°} form a right-handed orthonormal tetrad at 
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the observation point. Using the polarization basis vectors, the observer can decompose an 
arbitrary radiation field into Stokes parameters (e.g. Ref. /(i?,e°), Q{E,e°')^ U{E,e°') 
and V{E, e") along the direction at photon energy (or equivalently frequency) E. The 
transformation laws of the Stokes parameters under rotations of (ei)° and (62)" lead one to 
introduce a second-rank transverse polarization tensor Pab{E,e'^). The only non- vanishing 
tetrad components of Pab{E, e^) are 

for i and j = 1,2, and we have left the arguments E and implicit. Introducing the 
projected alternating tensor eabc = Vabcd.u'^, where rjabcd is the spacetime alternating tensor, 
we can write Pab{E,e'^) in the covariant irreducible form 

Pab{E, e'^) = -h{E, e^)nab + Vab{E, e'^) + W{E, e'^)e,fc,e^ (2.5) 
which defines the linear polarization tensor Vab{E,e'^), where 

'Pab{e.ne,f = \i^^ (2.6) 

which is transverse and trace-free. For a general second-rank tensor Sab, we follow 
Thorne |jl8| by denoting the transverse trace- free (TT) part by [S'ab]'^'^, so that 

[tiab\ = narif^Ocd — -rLabrL Ocd- {^■i) 

The magnitude (squared) of the linear polarization is + f/^ which can be written in the 
manifestly basis-independent form 2VabV°'^. 

Since I{E, e^) and V{E, e'^) are scalar functions on the sphere at a point in spacetime, 
their local angular dependence can be handled by an expansion in PSTF tensor-valued 
multipoles 

00 

I{E,e'^) = Y.lAAE)e^\ (2.8) 

(=0 

00 

V{E,e^) = Y.VAAEV^- (2.9) 

The expansion in PSTF multipoles is the coordinate-free version of the familiar spherical 
harmonic expansion. Eqs. (|2.8| ) and (O) can be inverted using the orthogonality of the 
g(^!>^ for example, 

lAXE) = Ar' JdQIiE,e'^)e^A,) where A, ^ (2.10) 

For the TT tensor Vab{E,e'^), we use the TT elements of the PSTF representation of 
the pure-spin tensor spherical harmonics (see Ref. |]18[ for a thorough review of the various 
representations of the tensor spherical harmonics), so that 



4 



1=2 1=2 



The first summation in Eq. ( p.ll[ ) defines the electric part of the linear polarization, while 



the second summation defines the magnetic part. (The electric and magnetic parts of the 
polarization were denoted by G(radient) and C(url) in Ref. [|^.) The l-th term in the 
multipole expansion of the electric part has parity (—1)', while the l-th term in the magnetic 
part has parity (— l)'"*"^. Eq. (|2.11|) can be inverted to give 

£A,iE)=Mi'Ar' I dne^A,.,Va,_,a,)iE,e'^), (2.12) 
-BaXE) = MMr'J dne,e'\a,eA,^,Va,_,)d{E,e^), (2.13) 



where M; = J2/(/ - l)/[(/ + 1)(/ + 2)]. The muhipole expansion in Eq. (|2]TT|) is the 



coordinate-free version of the tensor spherical harmonic expansion introduced to the analysis 
of CMB polarization in Ref. [jl9|, and the expansion in the Newman- Penrose spin- weight 2 
spherical harmonics employed in Ref. pO . 



A. Power spectra 

To define the bolometric power spectrum it is convenient to define energy-integrated 
multipoles, for example [| 

POO 

Ia,=Ai dElA^E), (2.14) 
Jo 



with equivalent definitions for Sai, Bai and Vaj- The factor A; is included in Eq. (|2.14|) so 
that the three lowest intensity multipoles give the radiation energy density, energy fiux and 
anisotropic stress respectively: 

J = p(^), Ia = q^^\ 4. = vri?. (2.15) 

In the almost-FRW case, with an ensemble that is statistically isotropic, the power spectrum 
of the (bolometric) temperature anisotropies, C/^, can be defined in terms of the covariance 
of the intensity multipoles 0: 

^,^2 



J J {w^n = ^icr6ih\2i, (2.16) 

where = hfa-^ ■ ■ ■ l^ai)- -^^^ small anisotropies, the gauge-invariant fractional temperature 
difference from the all-sky mean is related to the I a, by 



^Iai was denoted A^n^, in Ref. Q, and jj^^ in Ref. p. 
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-' 1=1 

so that the temperature correlation function evaluates to 

oo 

{ST{e'^)ST{e'^)) = E 



(2.17) 



1=1 



where X = — e'^e^, and Pi{X) is a Legendre polynomial. In deriving Eq. 



result e^^^'^e' 



(2.18) 



we used the 



(21 + l)A,PK^)/(4vr) 



The power spectra for the polarization are defined by analogy with Eq. (|2.16|) . We choose 
our conventions for the power spectra to conform with Ref. so that it is necessary to 
include an additional factor of M//\/2 for each factor of the polarization. For example, for 
the electric polarization 



') 



(/ + !)(/ + 2) 



(2.19) 



Unfortunately, the definitions of the polarization power spectra given in Ref. |2^ differ from 
those in Ref. |T9[ by factors of \/2; see Ref. and Sec. |VI| for details. For a parity 
symmetric ensemble the magnetic polarization is uncorrelated with the electric polarization 
and the temperature anisotropies. In this case, we find that the generalisation of Eq. ( |2.18| ) 
to linear polarization is 



{Va,{e^)V''\e'^)) = E 



1=2 



(2/+i: 

4:71 



crpi-i{x)], 



(2.20) 



where Pab(e^) = J^^ dE Vab{E , e") . When the directions e" and e'" coincide, Eq. (p:20|) 
reduces to the ensemble average of the square of the degree of linear polarization (expressed 
as a dimensionless bolometric temperature). Note that the electric polarization quadrupole 
gives an isotropic contribution to the linear polarization correlation function. Integrating 
Eq. ( [^.20D over all directions e''^, with e*^ fixed, the right-hand side reduces to 5Cf^. The 
presence of P/_2(X) and Pi^i{X) in Eq. ( p.20|) reflects the opposite parities of the electric 
and magnetic contributions to Vab{^^) for a given multipole /. 



B. Transformation laws under changes of frame 

The polarization tensor Pab{E, e^) has been expressed in 1+3 covariant form, so it is not 
invariant under changes of frame. If we consider a new velocity field = 7(u'^ + f"), where 
v"" is the projected relative velocity in the u"" frame and 7 is the associated Lorentz factor, 
for a given photon with 4-momentum p"" the energy and propagation directions in the 
frame are given by the Doppler and aberration formulae: 

E = 7E(l + eX), (2.21) 
~e- = [7(1 + e\)]-\u'' + e») - 7(m" + t;"). (2.22) 
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Note that e'^ is a projected vector relative to m". The screen projection tensor for a given 
null direction transforms to 

nab = nab - '^Pian.yv' + l-PaPbn^V%\ (2.23) 

while the polarization tensor transforms according to 

E-'Pab{E,e'') = E~'ntnt'Pd^d,{E,e''). (2.24) 

Under this transformation law the intensity, circular polarization and linear polarization do 
not mix. The irreducible components of Eq. ( p. 24 ) give the frame invariance of I{E, e'^)/ E'^ 



and V {E , e'^) / E'^ , and show that the transformation law for VabiE^e'^) follows that for 
Pab{e^)- The degree of linear polarization [2Vab{E^e'^)V°'^{E^e'^)Y/'^ / I{E^e'^) is invariant 
under changes of frame. The transformation law for Pab{E, e^) ensures that the tetrad com- 
ponents of Pab{E , e'^) / E^ are invariant if the polarization basis vectors are transformed as 

{e,r = ni{e.)\ (2.25) 

for 2 = 1,2. 

Under changes of frame, multipoles with different / mix because of Doppler beaming 
effects. To first-order in the relative velocity f the energy- integrated multipoles transform 
as 

4 = hf^)lB, -{I- 2)v%^, - (2.26) 

with an equivalent result for Va;. Although the linear polarization transforms irreducibly, 
the electric and magnetic parts do mix. To first-order in v°', we find 

f -UB,)^ m + 3) (/- 2)(/ -!)(/ + 3) , 
^A, - hi^Ai)^B, - (2/ + 1) (/Tip ' 

- ^^^Vbe'\a,BA,_,)c. (2.27) 



/(/ + 3) ^. ^ (/-2)(/-l)(/ + 3) , 

+ Ji^f^''\-^^^^^^> (2-28) 



In an almost-FRW model the polarization is a first-order quantity, as are physically defined 
relative velocities. It follows that the electric and magnetic multipoles are frame-invariant 
in linear theory. 

III. BOLTZMANN EQUATION 

The dynamics of the radiation field is described by the (exact) coUisional Boltzmann 
equation 
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C[E-^PabiE,e')] = Kab{E,e''), 

where the Liouville operator C acts on transverse tensors Aab = [Aab 
path x"(A), p"(A) in phase space, with = dx"/dA, according to 

C[Aab{E,e')] = nt'Kp'VeA,^,,{E,e''), 



(3.1) 

along the photon 
(3.2) 



where Va is the spacetime covariant derivative. The transverse scattering tensor Kab{E,e'^) 
describes interactions of the radiation with matter. The Liouville operator C preserves the 
irreducible decomposition of the polarization tensor [Eq. (|2.5| )], so that 



1 d 



C[E-'P^k{E,e'')] = -- — [E-'l{E,e'')]Hab + C[E-'Vab{E,e'')] 



2dA 
2dA 



[E-'V{E,e')]eabde' 



(3.3) 



In the absence of scattering, Kab{E,e'^) = 0, the tetrad components of Pab{E,e^)/E^ are 
constant along the photon path provided that the polarization basis vectors are transported 
along the null geodesic according to 



HabP^V^ie^f = 0, 



(3.4) 



and the constraint (cj) 

It follows from Eqs. (|2.24]) and ( |3.2| ) that under changes of frame H-i> u", the action of 
the Liouville operator transforms as 

C[E~^Pab{E,~e')] = nt"HtC[E-^Pd,d,{E,e'). (3.5) 

With this result, we deduce from Eq. ([Ol) that the scattering tensor Kab{E, e^) must trans- 
form as 



Kab{Ere') =K'HtK,,,,{E,e'). 



(3.6) 



This result is useful since the scattering tensor is often simplest to evaluate in some preferred 
frame, picked out by the physics of the scattering process. The scattering tensor in a general 
frame can then be computed using Eq. (|3.(j|). 

Over the epochs relevant to the formation of anisotropies and polarization in the CMB 
the dominant coupling between the radiation and the matter is Compton scattering. To 
an excellent approximation we can ignore the effects of Pauli blocking, induced scattering, 
and electron recoil in the rest frame of the scattering electron, so that the scattering may 
be approximated by classical Thomson scattering in the electron rest frame. Taking the 
electron 4-velocity to be m", and the (proper) number density of free electrons to be rze, the 
exact scattering tensor in the m° frame in the classical Thomson limit is 



E''kab{E,e')=h,aT 



\nab 



lU 5 



-^..(^,e-) + ^[4.(^)]^" + ^[4.(^)]^^' 



(3.7) 
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where ctt is the Thomson cross section. This expression for the scattering tensor follows 
from inserting the multipole decomposition of the polarization tensor into the Kernel for 
Thomson in-scattering (e.g. Ref. ^7^), and integrating over scattering directions. 

We have written Eq. (|3.7|) in irreducible form. The first set of terms in square brackets 
affects the evolution of the intensity in the electron rest frame, the second set of terms affect 
the linear polarization, and the third set affect the circular polarization. Concentrating first 
on the terms that affect the intensity, the lack of a monopole component is due to there 
being no energy transfer from Thomson scattering in the rest frame of the electron, while 
the presence of the intensity and electric polarization quadrupoles is due to the anisotropy 
and polarization dependence of Thomson scattering. Turning to the terms in Kab{E,e'^) 
which couple to the evolution of the linear polarization, the presence of the quadrupoles of 
the intensity and the electric polarization arises from Thomson in-scattering. These terms 
directly affect the evolution of the electric quadrupole alone in the electron rest frame, but 
this is not true in a general frame. Finally, we see from the terms in the final square bracket 
in Eq. ( |3.7| ) that the evolution of the circular polarization decouples from the intensity and 
the linear polarization. Since the circular polarization transforms irreducibly under changes 
of frame, circular polarization is not sourced by Thomson scattering for arbitrary frame 
choices. 



A. Integral solution for Vab{e^) 

The transformation of Kab{E,e'^) to a general frame gives a non-local expression in the 
energy E, since there is energy transfer in Thomson scattering from a moving electron. 
(See Ref. ||2^ for applications of this transformation in cluster physics.) It is therefore 
convenient to integrate over energy, in which case we find the following exact expressions 
for the evolution along the line of sight: 



POO A 

dEE^ — [E-'l{E,e'')] = -n,aT7(l + ^de'^)/(e^ 



+ -^neaT[7(l +t^de')]-=^(/ + aeV), (3.8) 
47r 



poo 

/ dE E^CiE-^'VabiE, 6^=)] = -neaT7(l + Vde'')Pabie' 
Jo 



+ l^n^^Tbil + vae')]-'[n:mi-Q,,,r^, (3.9) 



POO A 

dEE^ — [E-''V{E,e')] = -h^a^^^l + v,e'')V {e') 



- ^neaT[7(l + Vde'')]-^Vae\ (3.10) 
ovr 

The tildes denote that the quantity is evaluated in the rest frame of the electrons, which 

has four-velocity u"" = 7(ti" + f"), and we have defined e.g. I{e'^) = dE I {E , e'^) . The 
quantity (at, where 

3 9 

Cafe = -lab + -^ab, (3-11) 
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appears in Eqs. ( p.8D and (p.9|), evaluated in the u'^ frame, because of the anisotropy and 
polarization dependence of Thomson scattering. 

Formal solutions to Eqs. ( |3.8| - |3.10| ) can be obtained by integrating along the null cone. 
For the linear polarization we must work with the tetrad components of Pafe(e^), introduced 
in Sec. 0, since we can only meaningfully integrate scalar equations. Making use of 

{e,r{e,f / dEE'C[E-'VatiE,e')] = (e'^ + w")VjP,,(e^)(ei)'^(e,)^] 

^ 

+ {e^ne,fl dEE' — — [E-'Vat{E,e% (3.12) 

and integrating the second term on the right by parts, we find the following exact solution 
for the energy-integrated linear polarization tensor at some point i? in a general cosmological 
model: 

(3.13) 

Here, the integral is along the photon geodesic and the measure is dt = dx"'Ua. The redshift 
back from R along the line of sight is z, and f is the optical depth along the line of sight, 
measured in the electron rest frame: 

t\a = / dtnedT, (3.14) 

J A 

where dt = Uadx'^, and A is a point on the null geodesic through R. The integral solution, 
Eq. (|3.13|) , shows clearly how the intensity and electric quadrupoles generate polarization on 
scattering, and the effect of subsequent gravitationally-induced rotation of the basis vectors 
(cj)" as the radiation free streams to us. The need to calculate the intensity and electric 
quadrupoles, the redshift, and the rotation of the basis vectors along the line of sight suggests 
that in a general (non-linear) model it will be more convenient to work directly with the 
polarization multipoles. However, in Sec. |V] we show that the linearised version of Eq. ( p.l3| ) 
does provide a very direct route through to the integral solutions for the mode-expanded 
multipoles in the almost-FRW case. 



B. Multipole decomposition of the Boltzmann equation 



The Boltzmann equation ( p.l|) can be written in multipole form by expressing Pab{E, e'^) 
as a multipole expansion using Eqs. ( |2.8D , (|2.9|) , and (|2.11|) , and decomposing the resulting 
equation into multipoles. This leads to four sets of multipole hierarchies for IaX^), Sa^E), 
Bai{E), and Vai{E). In the absence of scattering the hierarchies for Ia^E) and V^, (£") 
decouple, but if Thomson scattering is included the hierarchy for Iai (E) includes the multi- 
poles of the linear polarization as source terms. The exact multipole equations for Iai (E) in 
the absence of scattering were given in Refs. ||7|, p!T| , p!2| ; the same equations apply to the mul- 



tipoles Vai of the circular polarization when there is no scattering. (The leading non-linear 
Thomson source terms in the hierarchy for the energy-integrated multipoles were also 
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given in Ref. 0, but these did not include polarization effects.) The exact multipole equa- 
tions for Sai{E) and Ba_,{E) are significantly more involved than those for the intensity or 
circular polarization. We defer the derivation of the exact polarization multipole equations 
to a subsequent paper. Here, we shall only consider the case of almost-FRW universes, in 
which case the polarization and I > intensity multipoles are 0(e) in a smallness parameter, 
e, which characterises the departure of the cosmological model from exact FRW symmetry. 
For the rest of this paper we shall restrict attention to almost-FRW models, and work only 
to first-order in e (linear perturbation theory). 

To first-order in e, we find for the electric polarization^ 

^ cutIBa,{E) = -n,aT£A,{E) + ^n,aT[Ia,a,iE) + 6£a,a,{E)]5f , (3.15) 



(/ + 1) "'^^ "'^^ 10 

and for the magnetic polarization 

BaSE) - \qe'^[E-'BaAE)] + liia,BA,_,){E) - lll]l[^lf B^AXE) 

^ cm\EAXE) = -n,a^BAXE)- (3.16) 



(^ + 1) 



In these equations, an overdot denotes the action of w^Va, is the totally projected 
derivative: 

T)aS^..., = Khl...h{V^S,...^, (3.17) 

for some arbitrary tensor 5'fc...c, and the volume expansion O = VqM". We have replaced 
by the electron number density in the frame, which is correct to the required order. The 
electric and magnetic multipoles are coupled through curl terms, where we have introduced 
the curl of a rank-/ PSTF tensor S'a,, 

curl^^, =e,,(,^D''^^^_^)^ (3.18) 

The coupling of electric and magnetic multipoles through curl terms is reminiscent of 
Maxwell's equations, as is the fact that only the electric multipoles have inhomogeneous 
source terms on the right. It is this curl coupling and the structure of the source terms that 
lead to the well-known result that scalar perturbations do not generate magnetic polariza- 
tion (see Sec. This is critical for the detection of a gravitational wave component to 



the anisotropy [P^]. 

It is often more convenient to work with the energy- integrated multipoles Eai and S^,, 
since it follows from Eqs. ( |3.15D and (3.16) that the polarization has the same energy spec- 



trum as the intensity anisotropies. Integrating over energies and then by parts, we find 



^We adopt the convention that Eai and Bai vanish for Z < 2, and Iai and Va, vanish for Z < 0. 
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c ,4 , (^ + 3)(/-l) , I 

' 3 ' — (l + 1)2 — ~ (2/ ' 



~ (7TT) " -^cC^T (^^A, - ^Caiaa*^/^) , (3-19) 

2 

+ curlf:^, = -neCrT^A;. (3.20) 

For the circular polarization, we find to 0(e) 

Va, + \qVa, - J^^^^iay^^,) + DV,^, = -ne^T (v^A, " ^K,^,^) , (3.21) 
and for the intensity 

= -neO-T (/a, - - ^/^^a,5/ - ^Ca,a,6^) , (3.22) 



which extends the result in Refs. [§,0,0 to include the polarization dependence of Thomson 
scattering. In Eq. ( p.22| ), Aa = Ua is the acceleration of and = D(a^6) is the shear. 
Eqs. ( p.l9| - pT^ ) provide a complete description of radiative transfer in almost-FRW models. 
They describe the evolution of the energy-integrated intensity and polarization multipoles 
along the integral curves of The equations are valid for any type of perturbation (we 
have not decomposed the variables into their scalar, vector or tensor parts), and for any 
(physical) choice of m". To close the equations it is necessary to supplement them with 
the 1+3 covariant hydrodynamic equations, e.g. Ref. |T3], to determine the 0(e) kinematic 
variables Aa, Va, and aab- 



IV. SCALAR PERTURBATIONS 

Up to this point our discussion has been quite general. Although we specialised to 
almost-FRW models when discussing the multipole propagation equations, we did not split 
the perturbations into their constituent modes (scalar, vector, tensor etc.). However, for 
detailed calculation with the linearised, almost-FRW equations, it is convenient to exploit 
the linearity to break the problem into smaller pieces. Our strategy follows Refs. P,P,pD|: 
we expand all 0(e) variables in PSTF tensors derived from appropriate irreducible eigen- 
functions of the comoving Laplacian S'^D'^D^, where S is the covariantly-defined scale factor 
with S/S = 6/3 and DaS = 0(e). Since the sources of anisotropy and polarization are at 
most rank-2 objects, there are non- vanishing contributions to the source terms only from the 
scalar, vector and rank-2 tensor eigenfunctions. The different perturbation types decouple 
at linear order, with each giving rise to a set of coupled, first-order ordinary differential 
equations. We consider the scalar modes in this section; tensor modes are treated in Sec. |V]. 
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Since vorticity dies away in expanding models in the absence of significant momentum densi- 
ties and anisotropic stresses we do not consider vector modes. If required, e.g. in defect 



models vector modes can be easily included in the formalism. The mode-expanded 1-1-3 
covariant equations describing the perturbations in the matter components other than the 
radiation, and the geometry, can be found in Refs. P,P,p!0[] , so we need only consider the 
radiation here. Furthermore, we do not consider circular polarization any further since it is 
not generated by Thomson scattering. 

For scalar perturbations we expand in the scalar eigenfunctions Q^^\ where 

S'^D"D„Q('=) = k'^Q^^^ (4.1) 

at zero order. The Q*^^-* are constructed so that Q^^"* = 0(e). The allowed eigenvalues /c^ 
depend on the spatial geometry of the background model. Defining z/^ = (fc^ -|- K)/\K\, 
where GK/S"^ is the curvature scalar of the spatial sections in the background model, the 
regular, normalisable eigenfunctions have u > for open and flat models [K < 0). In closed 
models u is restricted to integer values > 1 P^ , P7[] . The mode with u = 1 cannot be used 
to construct perturbations (its projected gradient vanishes globally), while the mode with 
z/ = 2 (which can only represent isocurvature perturbations |2^) only contributes to the 
CMB dipole. An explicit PSTF representation of the scalar harmonics Q^'^^ is given in the 
appendix. 

For the l-th multipoles of the radiation anisotropy and polarization we expand in rank-Z 
PSTF tensors, Qa\ derived from the scalar harmonics via 



The recursion relation for the Qa^, 

= |d(.,<_,), (4.3) 

follows directly from Eq. ([4.2|) . The factor of {S/kY in the definition of the Q^j ensures that 
Qa^ = at zero-order. For and Sai we writeQ 

/A.=/E«rVl')Qj/, />1, (4.4) 

k 

SA, = iT.^r'sPQAl (4.5) 

k 

where we have defined ai = ni=i f^n, with 

Ki = [l- - l)K/eYl\ I > 1, (4.6) 



^The ll ' are related to the J^ ' of Ref. [| by ij!' = aiJl >. 
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and tto = 1- The symbolic summation in Eqs. ( [4.4|) and ( [4.5|) denotes a sum over the 
harmonics (see Eq. ||A71| ), and the mode coefficients, such 0(e) scalars with O(e^) 

projected gradients. We need not consider the magnetic polarization since it is not generated 
by scalar modes [p!9| , p0| . To see this, we need the ffist-order identity 



curlD(„,^^^^,) 



(^ + 1) 



D(a,+i curl^Aj, 



(4.7) 



where Sai is an 0(e), rank-/ PSTF tensor. Applying this identity repeatedly to the right- 
hand side of Eq. (|4.3|) we see that the curl of Q^^^ vanishes, as we would expect for a 
scalar perturbation. It follows that in linear theory curl£^A, = for scalar perturbations, 
so the inhomogeneous source terms in the propagation equation for Bai (Eq. ||3.2CI|| ) vanish. 
Since primordial polarization is erased by Thomson scattering during tight coupling, scalar 
perturbations do not support magnetic polarization. 

The decomposition of I{e'^) into angular multipoles 7^,, and the subsequent expansion 
in the Qa\ combine to give a normal mode expansion which is equivalent to the Legendre 



tensor approach, ffist introduced by Wilson [|15|. The advantage of handling the angular and 
scalar harmonic decompositions separately is that the former can be applied quite generally 
for an arbitrary cosmological model. Furthermore, extending the normal mode expansions 
to the polarization and tensor modes is then trivial once the angular decomposition has been 
performed. 

Substituting Eq. (U) into Eq. (|3]2|), and using the zero-order identity 



k I 

S{21-1) 



Q 



(4.8) 



gives the Boltzmann hierarchy for the 



(0. 



k 

s 



{21 + 1] 



{21 + 1) 



til I, 



(i-i) 



+ 



4 ^ . r-1 8 k ^2 / t(i) 

+ -^g^kSi - —-K2ak6i = -nedT ( I'k 



It should be noted that ij^^ refers to the projected gradient of intensity monopole /, 

k 



3 



3^' 



(4.9) 



(4.10) 



which is gauge- invariant. To obtain Eq. ( [4.9| ) for / = 0, take the projected gradient of 
Eq. ( p. 22 ) for Z = and commute the derivatives. This gives rise to the source term Zk, 
where 



(4.11) 



The other kinematic source terms come from the acceleration Aa, the the baryon relative 
velocity Va, and the shear aab, with 
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^VkQ 



(fc) 

a 1 



k 



(k) 



(4.12) 
(4.13) 

(4.14) 



The term describing the dependence of Thomson scattering on the anisotropy and polariza- 
tion is Ck = ^Ik^/4: + 9£k^/2. Eq. (|;9|) extends the 1+3 covariant results of Refs. [|,| by 
including this polarization source term. 

For the electric polarization we use Eq. ([4.5|) in Eq. (|3.19|) to find 



k 
s 



(/ + 3)(/ 



(2/ + l)(Z + l) 



l^l+l^k 



{21 + 1) 
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CkSf 



(4.15) 
and 



which is the 1+3 covariant analogue of equivalent results in Refs. |p!4| , p9| . Eqs 
( [4.15| ) can be integrated with 1+3 covariant hydrodynamic equations |^,|] to determine the 
anisotropy and polarization from scalar modes. 

To obtain formal integral solutions to the Boltzmann hierarchies we note that the homo- 
geneous form of Eq. ( [4.9| ), obtained by setting ne = and removing the kinematic source 



terms, is solved by $f (x), where x = y\K\{ri[i — rj), with rjn the conformal time at our cur- 
rent position R. Here ^'{{x) are the ultra-spherical Bessel functions [^. The full solution 
of Eq. (^4.9| ) follows from Green's method: 



r(0 



tn { ( k 1 ^ 



k 



vV2+Tdx 



^x) - 



1 fk 
3 



(z/2 + 1) da;2 



Zk - QAk 



<|.r(x) 



(0) 



$r(x) , (4.16) 



where r is the zero-order optical depth back to x (see Eq. ||3.14|] ). The geometric factors 
$['/3 + (z/2 + l)~M2$^/dx2 and (z/^ + l)-^/'^<\^'{ /<\x arise from the projections of QlJe'^e'' 
and Q^^e'^ respectively, at x back along the line of sight. 

For the polarization we note that the homogeneous part of Eq. (|4.15|) is solved by /(/ — 
l)$^(x)/ sinh^x. The full solution is 



(I) 



(z/2 + 1) 



smh X 



(4.17) 



Eqs. ( [4.16|) and ( [4.17| ) are valid in an open universe. For closed models one should replace 
the hyperbolic functions by their trigonometric counterparts, and z/^ + n by z/^ — n where n 
is an integer. 

To characterise the initial amplitude of the mode we introduce a set of random variables 
(pk with the covariance matrix 



) 



|i^|3/2 ^(j,2 + 1) 



Jkk' 



(4.18) 
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appropriate to a statistically isotropic and homogeneous ensemble. The 5kk' is defined by 
J^k'^k^kk' = Afc/ for any mode functions (see appendix). We follow the conventions of 
Lyth & Woszczyna ^0[, so that the scalar field Y,k 'PkQ^'^^ has a scale-invariant spectrum for 
P^(i/) = constant. We write the radiation mode coefficients in terms of the (p^ and transfer 
functions Tf\v) and T^\v): 

ir = TPw)A. 4" = ^ri'V)^», (4.19) 

where the factor M;/a/2 is for later convenience. To compute the anisotropy and polarization 
power spectra we use Eqs. ( [4.4| ) and ( [4.5|) in Eqs. (p.l6|) and ( p. 191) to find 



Cr = i,L 7:^,T^\-)T?iu)V,{u). (4.20) 



1 udu 
leio (^ 

Here, X and Y represent either of / or In a closed universe the integral over u in Eq. ([4.20 ) 
should be replaced by a sum over integral u > I + 1. To derive Eq. ([4.20|) we have made use 
of the result 

E f{r^)Q'iy'''^' = J^MiApi [ \Kf-v^f(v)a,\ (4.21) 
for any scalar function /(z^), which follows from Eqs. (|A6|) and (|A9|) in the appendix. 



V. TENSOR PERTURBATIONS 



For tensor perturbations we follow the procedure in Ref. [|T0| and expand in rank-2 PSTF 



tensor eigenfunctions Qit'' of the comoving Laplacian: 



(5.1) 



The tensor harmonics are transverse, D'^Q^^^ = 0, and constant along the integral curves of 
m", Qq^^ = 0, at zero-order. For tensor modes we define = {k'^ + 3K)/\K\. The regular, 
normalisable eigenmodes have u > for fiat and open models, while for closed models u is an 
integer > 3. Explicit forms for the tensor harmonics in the PSTF representation are given 

ITUi 



in the appendix; see also Ref. []TU| . The tensor harmonics can be classified as having electric 
or magnetic parity. We continue to denote the electric parity harmonics by Q^^b ^ but we use 
an overbar to distinguish the magnetic parity harmonics: Q^ab ■ The electric and magnetic 
parity harmonics are related through the curl operation (see Eqs. [|A19|| and ||A20|| ). 

Following our treatment of scalar perturbations, we form rank-/ PSTF tensors Q^^^ and 

^ (k) 

Q]^ from the electric and magnetic parity tensor harmonics respectively. We define 



/-2 



D(ai • • • Da,_2 



Q 



(k) 



(5.2) 
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with an equivalent definition for the magnetic parity harmonics. The rank-/ tensors satisfy 
QV = at zero-order. We expand the /-th multipoles of the intensity and hnear polarization 



as 





k 


Vk ^Ai + -'fc VA( J? 


I > 2, 


(5.3) 




k 


{<-'k ^Ai + <-'k ^Ai J' 




(5.4) 


^A, 








(5.5) 



For tensor modes we have defined Pi = ni=2 ^n, and 



211/2 



/ > 2. 



(5.6) 



It is necessary to include the magnetic polarization for tensor modes since the source to the 
magnetic hierarchy does not vanish: cml Sa, 7^ 0. 

The Q^^l constructed from the tensor harmonics satisfy the same recursion relation, 



Eq. (|4.3| ), as their scalar counterparts. However, their projected divergences evaluate to ||T0 



5/(2/- 1) 



l-(^'-3)- 



A;2 



(5.7) 



An equivalent relation holds for the Qa^- Using these results in Eq. ( |3.22| ) we find 



r(0 



k 

s 



(/ + 3)(/-l) 



/ 



" (27TT)"'^' 



(/-I) 



8 k 
155 



,(/ + l)(2/ + l) 
2 

15' 

where (k = S/^^V^ + Qi^^^V^, and cr^ is the tensor mode coefficient for the shear. 



-neCTT ( - -^CkSf I , 



(Tab = J2 ^i^kQah + ^^kQab)- 
k 



(5.8) 



(5.9) 



Eq. (|5.8|) extends the 1+3 covariant result in Ref. [|T^] to include the polarization source 
term. For the polarization hierarchies we need the result 



curlQV = y-Wl + — gy. 



which follows from repeated application of Eq. ( ^4.7| ), and finally Eq. 
curlQ^^^^^^ in terms of Q^J^^-^ar For the electric polarization, Eq. ( p.l9| ) gives 



(5.10) 
to express 



s 



(2/ + l)(/ + l)3 



(2/ + 1 



3K 



l{l + l)S 



1 + ^^. 



(0 



(5.11) 
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Eq. ( p.20|) gives the corresponding equation for the magnetic polarization: 



S 



(/ + 3)2(/ 



(2/ + l)(Z + l) 



k 



{21 + 1) 




l{l + l)S 



(5.12) 



Equivalent equations hold for the barred variables 8^^ and B'^^' . Note how power is trans- 
ferred between B^l^ and 8^^ due to the curl coupling between the electric and magnetic 
multipoles in Eqs. ( ^.191 ) and ( |3.2CI| ). 

The integral solution for the tensor contribution to the intensity anisotropy in a general 
almost-FRW model was given in Ref. Including polarization modifies the scattering 

source term, so that the integral solution becomes 



?(0 



r(0 



4/(Z - 1) 



[(z^2 + l)(zy2 + 3)]l/2 



k 

r 



1 

-1 
4 



sinh^x' 



(5.13) 



in an open universe. The geometric factor /(/ — l)$^(a;)/ sinh^x follows from the projection 

(k) h 

of Qaf, e°e at x back along the line of sight. 

To solve the coupled equations (|5.11|) and (|5.12| ) for the linear polarization it is simplest 
to return to the integral solution for the tetrad components of Vab{e'^), given as Eq. ( p. 131) . 



In linearised form, we have 



(5.14) 



Substituting the tensor mode expansion of C,ab gives terms in the integral like 
[Q^ahY'^ i^iYi^iY from the electric parity harmonics, and similar terms from the harmon- 
ics with magnetic parity. Using the representation of the tensor harmonics given in the 
appendix, we find that 



[Q 



T3(x)[Q, 



(5.15) 



where T^ix) is a z/ and L-dependent function given by Eq. (|A13| ). The superscript [k) on the 
tensor harmonics represents the collection z/, L, and M, where L and M describe the orbital 
angular momentum of the harmonic; see the appendix for details. The Q^^^^^ are rank-L 
PSTF tensor fields, which are introduced in the appendix. They encode the 2L + 1 degrees 
of freedom in the tensor harmonics with L units of orbital angular momentum. Using the 
transport properties of Q^^^^\ Eq. ( |A3D , and of the tetrad vectors (cj)", Eq. (p.4|) , it is 
straightforward to show that to zero-order 



(5.16) 



so that [QlS_,e^-^]^^(e. 
harmonics, the representation in the appendix gives 



{ejY is constant along the line of sight. For the magnetic parity 



[QW\e,ne,f = f,{x)[e^,. 



(i^b)Ci_i^ 



iTT, 



(5.17) 
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where T'i(x) is given by Eq. ( |A18|) . Working to zero-order, we find that 

{u^ + e")V,{[e,,e^-'=--(,Qi[*J^^e^--]^^(e.)'^(e,)n = 0, (5.18) 

so that [gcl^'^^^'^'^ {aQb)c^j^_-^^^^~^V^ i^iYi^jY is also constant along the line of sight at zero- 
order. With these results, Eq. ( p. 141) reduces to 

k 

+^Ebc.e^-^-^(,Qi|*^)^e^-^]^^/*"dtneaTe--«:2-^aTi(x). (5.19) 

Extracting the Z-th electric and magnetic parity multipoles from this equation, we find 

^a\r = ^E^'i^'2?f V f''dtn,aTe-^K2-'aTs{x), (5.20) 



k 

^a,\r = ^T.^ilAiQ^a''''\r f''dtn,aTe-^K2-'(kf,{x). (5.21) 

To determine the harmonic coefficients, such as sl^\ at we note that Q^^^ = at i?, so that 
only the electric parity harmonics contribute there. We can use Eq. ( |A15D to substitute for 
SilQa'I^'' each u in terms of QaiIr in Eqs. ( |5.20| ) and ( |5.21| ). Comparing with Eqs. ( p.4| ) 
and ( |5.5| ), we can read off the integral solutions 

4') = Ml r dtn,aTe-^K2-'CkT3{x), (5.22) 

en 1 



Bt> = Mi^== dtn,aTe-^K2-XkTiix). (5.23) 
Evaluating the functions T^lx) and Ti(x), we can write the integral solutions in the form 



M,2 ft 



^^'^ " 2[(^2 + l)(!/2 + 3)]i/2 j dtn^aTC 'Ck(pi{x), (5-24) 



2[(z/2 + l)(z/2 + 3)]l/2 

where the geometric terms 0^(x) and ipii^) ^-re 

d2 d 
0J'(x) = — $;'(x) +4cothx— $r(x) - (z/2 - 1 - 2coth2x)$r(x), (5.26) 

(Jiilj Lit// 



^r(a;) = -2z/ 



-5^<l>r(x) + 2cotha;$r(x) 
dx 



(5.27) 



in an open universe. In closed models we replace i/^ + n by z/^ — n, and the hyperbolic 
functions by their trigonometric counterparts (as for scalar perturbations). It is straight- 
forward to verify that the integral solutions, Eqs. ( p.24| ) and ( |5.25|) , satisfy the multipole 
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and bIIi^ but with the replacement Cfc ^ Cfc for the source function. Although the magnetic 
parity tensor harmonics do not contribute at R, they do contribute to the anisotropy and 
polarization at points not on the integral curve of that passes through R. The integral 
solutions given here agree with those of Ref. \T^\ which were derived with the total angular 
momentum method. 

For tensor perturbations we characterise the initial amplitude of the modes by introduc- 
ing random variables (pk and (pk with the covariance structure 



equations and ( |5.12|) . The solutions for S^^'' and B^^^ are the same as those for S,^ 



) 



^k<Pk 



) 



1 



-6, 



kk' 



■)k(Pk 



) = o, 



(5.28) 



appropriate to statistical homogeneity and isotropy. Note that this form for the covariance 
structure forbids any cross correlations between the Eai and £ai or 7^,. In Ref. [|10| we 
took 0fc = (^^ + 3)£'/c/(z/^ + 1) where Ek is the initial amplitude of the dimensionless mode 
coefficient representing the electric part of the Weyl tensor. In this case, the minimal scale- 
invariant prediction of one bubble open inflation is = tanh(7rz//2) ||31]]. We write the 
mode coefficients for the intensity and the polarization in the form 



r(0 



Ti'\iy)<Pk 



4'^ = ^T«, 



( I ) _ Ml ( I ) - 

~ "71 



(5.29) 



Note that B^^:^ is proportional to (f)k since the source to B^!^ is (k rather than (k, see Eq. ( p.25p . 
We can now compute the anisotropy and polarization polarization power spectra for tensor 
modes from Eqs. (|2.16| ) and ( p. 191) . The result is 



cr 



1 (/ + 2)(/ + l) 
16 2/(/ - 1) 



z/dz/ (z/2 + 3) 
fz/2 + 1) 7^ 



(5.30) 



where XY is equal to //, S£, BB, or IS. In closed models the integral over u is replaced by 
a discrete sum over integral u > I + 1. In deriving Eq. ( |5.30|) we used the following result: 



k 



{k)^(k)Bir 



(47r) 



(5.31) 



for any scalar function f{v). This result is easily verified at R by using Eq. ([A15| ) and the 
fact that Q^S\r = 0. 



VI. DISCUSSION 



Most of the modern literature on the polarization of the CMB employs representations 
based on either a coordinate representation of the tensor spherical harmonics, introduced to 
CMB research in Ref. |T9|, or the Newman- Penrose spin- weight 2 harmonics (e.g. Ref. [0), 
first used for describing polarization in Ref. ||2^. Given the widespread use of these two 



formalisms, it is worthwhile outlining their relation to the PSTF representation adopted 
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here. We begin by introducing an orthonormal triad of projected vectors {(71)°, (72)", (73)°} 
at the observer's location. At this point, the triad is used to define angular coordinates 
{6, 0} which cover the two-dimensional manifold of unit projected vectors {n°}. The scalar 



spherical harmonics Y(^im){ 
3^i'r^' where 



Y{im){n°') can be represented by complex PSTF tensors 



A, 



dQY( 



(Im) 



n ]n 



(6.1) 



so that Y( 



(Ini) 
■.(Im) 



n 



yA^ 



n 



Using the spherical harmonic addition theorem, one can show 



that the y)^^ satisfy the orthogonality relations 



A 



(6.2) 
(6.3) 



which are useful for the discussion below. 

If we compare the PSTF expansion of the temperature anisotropy, Eq. ( p.l7| ), with 
Eq. (2) of Ref. |]19| we can read off the relation between the PSTF multipoles 1^, and the 
complex scalar multipole coefficients a^^^ of the anisotropy: 



TT 



lA, 



'-lYL 



(6.4) 



The factor of (—1)' in this equation arises because in Ref. the unit projected vector 
represents a direction on the sky, rather than the photon propagation direction, so that 
gO _ Using Eqs. ( |0|) and ( |6.4|) in the expression for the temperature power spectrum, 
Eq. ( pl6D , we find that 



/ T* T \ 



(6.5) 



It follows that the of Ref. |]T9[ coincides with the temperature power spectrum defined 
here. 

For the linear polarization we need to relate the PSTF representation of the TT tensor 
spherical harmonics to the coordinate representation derived by covariant differentiation 
of the y*^'™^(n"). Taking covariant derivatives on the sphere of the PSTF representation 
of the F('™)(n"), and using the conventions for the G(radient) and C(url) tensor spherical 
harmonics of Ref. 19 , we find that 



„Ci 



abC, 



n 



iTT 



,(im) 



Ci-2lTT 



(6.6) 
(6.7) 



Integrating Eq. (p.ll|) over energy and dividing by vr/J expresses the linear polarization in 
dimensionless temperature units. Comparing with Eq. (2) of Ref. |]T9|], we find the following 
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relations between the coordinate-dependent, scalar- valued multipole coefficients a^^) and 



defined there, and the PSTF multipoles and employed here: 



TT 

7a; 



m=—l 



TT 

7a; 



(6.8) 
(6.9) 



m=—l 



Using these results in the definitions of the polarization power spectra, e.g. Eq. ( p. 191) 
we find the non-vanishing covariance structure 



2{a 



c 



G* 



Cf^ Sill 5mm', 
Cf^5lli5mm'-i 
Cl^6ll'5mm' ■ 



(6.10) 
(6.11) 

(6.12) 



we find that 2Cp 



in a parity-symmetric ensemble. Comparing with Eq. (8) of Ref. [19 
and 2Cp = Cf^, while v^C^G = C/^. 
To compare the PSTF representation of the linear polarization with the spin-weight 2 
representation of Ref. 
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it is simplest to express the coordinate representation of the 
tensor spherical harmonics of Ref. |T9| in terms of the spin-weight 2 spherical harmonics. 
Noting that the Stokes parameters in Ref. pO| are expressed on a right-handed basis with the 



3-axis opposite to the direction of propagation (which fiips the sign of U), we find that for the 



ttE^im and as^im multipole coefficients of Ref. ||20|, a^;^ 



Im 



G 

(Zm) 



and as. 



Im 



whereas the temperature multipoles are equal. It follows that the polarization power spectra 
of Ref. pO[ are related to those defined here by Cei = Cp, Cbi = Cf^, and Cci = —Cj^. 



VII. CONCLUSION 

We have introduced a new multipole formalism for describing polarized radiation on the 
sky. In this approach the polarization tensor is expanded in the PSTF representation of the 
transverse-traceless tensor spherical harmonics . The PSTF representation is particularly 



convenient since the radiation multipoles are then coordinate-independent. This allows the 
equations of radiative transfer in a general spacetime to be recast as multipole equations, 
which is essential for the transfer problem in optically thin media. We applied the formalism 
to give a rigorous discussion of the generation and propagation of CMB polarization in 
cosmological models. We gave new results for the non-linear transformation properties of the 
polarization multipoles under changes of reference frame, the exact source term for Thomson 
scattering, and the multipole propagation equations in linearised form in an almost-FRW 
model. It was not necessary to split the perturbations into scalar, vector and tensor modes 
to obtain these results, and so they provide a solid foundation on which to build a complete, 
second-order analysis of CMB polarization. By expanding the linearised multipole equations 
in scalar and tensor harmonics we derived the mode-expanded multipole equations in almost- 
FRW models with general geometries, and their integral solutions. These results confirm 



those obtained earlier by Hu et al. with the total angular momentum method. 
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The linearised results of this paper have been included in a publically available Fortran 
90 code 1 3^, based on CMBFAST |1^, for calculating the CMB temperature anisotropy and 



polarization in general FRW models, including those with closed geometries, within the 1+3 
covariant and gauge-invariant approach. 

In a subsequent paper we will complete the development of a multipole transfer formal- 
ism in general spacetimes by providing the exact polarization multipole equations for an 
arbitrary geometry. This should simplify the modelling of a number of important astro- 
physical situations where relativistic flows or strong gravity effects are important, including 
non-linear effects in the CMB. 
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APPENDIX: SCALAR AND TENSOR HARMONICS 



In the text we make use of a specific PSTF representation of the scalar and tensor har- 
monics. This representation was described in detail in Ref . [|l^ , so will only be summarised 
here. The starting point is to "coordinatise" the solution in terms of a projected vector field 
e"", and a scalar field x- The restriction of these fields to the past lightcone through some 
point R (which can conveniently be taken to be our point of observation) are the propaga- 
tion direction of a free-streaming photon which will subsequently pass through R, and the 
conformal look-back time along the photon path respectively. The projected fields e"" and 
X are generated from their restriction to the past lightcone by Fermi transporting along the 
integral curves of u": 



e 



=0, X = 0. (Al) 



It is straightforward to show that the field satisfies e'^DaC^ = at zero-order in an almost- 
FRW universe. 



1. Scalar harmonics 

In the PSTF representation the regular, normalisable scalar harmonics can be written as 

Q(') = $I(a;)Q;,f ^e^^ L > 0, (A2) 

where the Q^^^^^ with M = —L . . . L are rank-L PSTF tensor fields satisfying the zero-order 
equations 

e'B.Ql''^ = 0, Qji^) = 0, (A3) 

which determine the fields from their initial values at R. The ^l{x) are ultra-spherical 
Bessel functions (see, e.g. Ref. [^) with z/^ = (/c^ + K)/\K\ and x = ^/\K\x■ In closed 
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models u is restricted to integer values with v > L. For convenience we label the harmonics 
with a lumped superscript {k) which represents z/, L and M. 
The scalar harmonics are normalised so that 

d^]e^l,dx sinh2a:g('=)g(^') = |5^^.A,Q!^f - u'), (A4) 

where di^eii^ denotes an integral over solid angles at R. Eq. (^) refers to an open universe. 
In closed models the hyperbolic functions replace their trigonometric counterparts, 5{v — u') 
replaces S^yi, and the upper limit on the normalisation integral becomes tt. It is very 
convenient to choose the Q^^^"* at R so that 

Q!ff^Q(^*^')^^ = A^-^W, (A5) 

which implies 

E Q!£^)Q(^*^)^-=A,-^/^g^J. (A6) 



M=-L 



It follows from Eq. (|A3|) that Eqs. ( [A5|) and ([A6|) hold at zero-order at all points. A 
scalar- valued, statistically homogeneous random field, say A(x"), can be constructed from 
a superposition of the scalar harmonics: 

A(a;'^) = / \K\'/'uMu ^ E ^uLmQ^'\ (AT) 

L=OM=-L 

in an open universe. In the closed case the integral is replaced by a discrete sum over integer 
I', and the sum over L is restricted to L < u. In the text we denote the sum over scalar 
harmonic modes in the symbolic form A(x") = Y.k^kQ^''^- The covariance structure 

(AfcAfe,) = A\u)Skk' (A8) 

is sufficient to ensure statistical homogeneity and isotropy of the ensemble. The symbolic 
Skk' represents 5ll'Smm'\K\^^^'^i'^^S{i' — u') in open models. For closed models S{u — v') 
should be replaced by 8^1,1. 

From the Q^^^ we can derive rank-/ PSTF tensors Q^^^ as in Eq. ( [4.2| ). For the specific 
representation given here we find that at the point i?, 

Qt\R = ^^mQ^A''^^iL. (A9) 

so that the only modes to contribute to the Z-th multipole of the radiation anisotropy and 
polarization at R have / units of orbital angular momentum. In a closed universe, the 
vanish globally for / > v. 
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2. Tensor harmonics 



a. Electric parity 

The regular, normalisable tensor harmonics with electric parity are given by 

+T,{x)e^^Hl^Q^cT^e''^-^ + T,{x)[Qii^'l^e'^--^f\ L>2. (AlO) 

Here, the (screen) projection tensor Hab = hat + e°'e'', and the subscript TT denotes the 
transverse (to e"), trace-free part. In an open universe, Ti{x) is given by 



Tiix) ^ 



{L + 2)\ 
(L-2)! 



1/2 



sinh^x ' 



(All) 



z/^2(z/2 + 1 
and T2{x) and T^{x) are determined by 

The electric parity tensor harmonics are normalised so that 

/"dfi,»l,dxsinh\Ql?Q('=')'^'' = ^5l^.A,.Q^/^)q(^*^')^-z.-25(z. - z.')- (A14) 

For closed universes, the hyperbolic functions should be replaced by their trigonometric 
counterparts, v'^ + n should be replaced by z/^ — n with n an integer, and 5{v — v') should 
be replaced by 5^^'. For closed models, the regular, normalisable modes have v an integer 
> 3, restricted to u > L. 

We can form rank-/ (/ > 2) PSTF tensors from the electric parity tensor harmonics 
as in Eq. ( ^.2|) . At the point R we find that 



qSIk = f ^^M-A.AQir^.., (A15) 

Air u 

so that, as with scalar modes, the only contribution to the l-th multipole of the anisotropy 
and polarization at R comes from those modes with / units of orbital angular momentum. 

(k) 

In closed models the vanish for I > u. 

h. Magnetic parity 

For the magnetic parity tensor harmonics, which we denote with an overbar, the regular 
normalisable solutions are 
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Ti{x)[e, 



q{LM) 



iTT 



(A16) 



with L>2. In an open universe, T2{x) is given by 



T2(X) 



(^+1) 



ni/2 r 



iL-2y. 



1/2 



sinhx 



and Ti{x) is determined by 



(A17) 



TAx) 



{L + 2) sinh X da; 



[sinh^a;T2(a;)]. 



(A18) 



In a closed universe we make the same replacements as for the electric parity harmonics. 
The magnetic parity harmonics satisfy the same normalisation condition, Eq. ( A14| ), as the 
electric parity harmonics, to which they are orthogonal. We can also form rank-/ PSTF 
tensors Q^^^ from the magnetic parity harmonics, as in Eq. ( ^.21) . However, we now find 

that the vanish at R, so that the magnetic parity harmonics do not contribute to the 
anisotropy and polarization there. (However, magnetic parity modes do contribute at points 
not on the integral curve of that passes through R.) In closed models the Q^^^ vanish 
globally for I > v. 

The electric and magnetic parity tensor harmonics are related through the curl operation 
(which is parity reversing). For the conventions adopted here, we find 



curlQi? 
curlQi? 



k 
S 

k I 3K 



1 + 



(A19) 
(A20) 



c. Statistically homogeneous tensor-valued random fields 

We can construct a statistically homogeneous and isotropic tensor-valued random field 
by superposing the electric and magnetic parity tensor harmonics. For example, for the 
shear tensor aab, we write 

L 



aab= TlKl'/'iy'duY: E ^{^uLmQ'S + CT.LMQ'ab^ (A21) 

L=2M- ^ ^ 



in an open universe. (Equation ||A21|| can easily be generalised to include supercurvature 
modes, if these are present in the initial conditions.) In a closed universe the integral over 
u is replaced by a sum over integer z/ > 3, and the sum over L > 2 is further restricted to 
L < u. In the text we use the symbolic summation J2k represent the sum over modes on 
the right-hand side of Eq. (|A21|) . The requirements of statistical homogeneity and isotropy 
restrict the covariance structure of the ak and (Xfc to the form 
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{akCTk') = 0, (A22) 

where the symbohc 6kk' is the same as for scalar modes. The shorthand ak is used to 
represent a^iM, and similarly for the barred variables, in Eq. ( |A22| ) and also in the text. 
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